The consequences of large #13 for the turbulence signatures in supernova neutrinos 
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The set of transition probabilities for a single neutrino emitted from a point source after passage 
through a turbulent supernova density profile have been found to be random variates drawn from 
parent distributions whose properties depend upon the stage of the explosion, the neutrino energy 
and mixing parameters, the observed channel, and the properties of the turbulence such as the 
amplitude C*. In this paper we examine the consequences of the recently measured mixing angle #13 
upon the neutrino flavor transformation in supernova when passing through turbulence. We find the 
measurements of a relatively large value of #13 means the neutrinos are relatively immune to small, 
C* < 1%, amplitude turbulence but as C* increases the turbulence effects grow rapidly and spread 
to all mixing channels. For C* > 10% the turbulence effects in the high (H) density resonance mixing 
channels are independent of #13 but non-resonant mixing channels are more sensitive to turbulence 
when 0i3 is large . 

PACS numbers: 47.27.-i,14.60.Pq,97.60.Bw 



I. INTRODUCTION 



The progress in the field of supernova neutrinos over 
the past decade has frenetic. The rich phenomenolo gy o f 
neutrino collective effects [i"14l5| (for a review see [l6lll7l|) 
have received the most attention but there has been 
an equally radical overhaul of the Mikheyev, Smirnov 
& Wolfenstein (MSW) [H, effect as applied to su- 
pernova ever since it was realized by Schirato & Fuller 
[2fj | that the Shockwave racing through the stellar man- 
tle could leave an imprint upon the neutrinos emitted 
from the cooling proto- neutron star (20rl26| . Recent stud- 
ies indicate Earth matter effects may be minimal [27| . 
From this ever-growing body of literature one now ex- 
pects that the neutrino signal from the next supernova 
in our Galaxy will be pregnant with information. If the 
signal can be decoded we might be able to both determine 
any unresolved properties of the neutrino and also to ob- 
serve the explosion while it is still deep within the star. 
Yet most, though not all, of these studies use spherically 
symmetric density profiles either in a parametrized form 
or taken from one-dimensional hydrodynamical simula- 
tions. While the use of one-dimensional hydrodynami- 
cal profiles for neutrino signal construction is probably 
adequate for certain situations - such as neutrinos from 
Oxygen, Neon, Magnesium supernova [2814301 which ex- 
plode in spherically symmetric simulations [3ll433| - it is 
now apparent that iron core collapse supernova should 
not be expected to be spherically symmetric. Large scale 
inhomogeneities are created deep within the explosion 
and one observes turbulence during the neutrino heat- 
ing/Standing Accretion Shock Instability phase [4ll448| 
leading to the expectation of violent fluid motions and 
turbulence in the mantle as the shock is revived and 
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moves outwards. Like collective and shock effects, tur- 
bulence is another supernova feature that can leave its 
fingerprints upon the neutrino burst [10, 49-5§|. At first 
glance turbulence is just a case of a more complicated 
MSW effect but, upon further reflection, one realizes that 
the randomness of the profiles means the transition prob- 
abilities for a particular neutrino - the set of probabili- 
ties that relates the initial state to the state after passing 
through the supernova - along a given ray are not unique: 
they will depend upon the exact turbulence pattern seen 
by the neutrino as it travelled through the supernova. 
The transition probabilities are drawn from a distribu- 
tion whose properties will depend upon the stage of the 
explosion, the character of the turbulence, and the neu- 
trino energy and mixing parameters. When the mixing 
angle #13 was unknown it was difficult to make robust 
statements about the effect of turbulence because at one 
value of #13 the effects would be negligible, at another 
the turbulence would be endemic. The recent measure- 
ments of the last mixing angle #13 by T2K [34| , Double 
Chooz H3, RENO [H and Daya Bay [H are all in the 
region of 9^ ~ 9°, significantly higher than the Dighe & 
Smirnov |38| threshold, and it is now possible to be more 
definitive about the consequences of turbulence. 



In this paper we consider the implications of the recent 
measurement of #13 upon the neutrino transition proba- 
bilities as a function of the turbulence amplitude. Our 
calculations expand upon the work of Kneller & Volpe 
[52| upon which we shall rely heavily for the techniques 
used to calculate the turbulence effects and as reference 
for our results. We first describe the calculations we un- 
dertook then present our results for the turbulence effects 
when the turbulence amplitude is small, less than 1% 
comparing large and small 6*13. We then turn to large am- 
plitude turbulence and compute the expectation values of 
the transition probability distributions in both neutrinos 
and antineutrinos again comparing large and small #13. 
We finish with a summary and our conclusions. 
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II. DESCRIPTION OF THE CALCULATIONS 

The quantities we are interested in calculating are the 
probabilities that some initial neutrino state \v(x)) at x 
is later detected as the state \v{x')) at x' . These prob- 
abilities are computed from the S-matrix which relates 
the initial and final states via \v(x')} = S(x',x) \is(x)). 
The S'-matrix is found by solving the equation 



where H is the Hamiltonian. In matter the Hamiltonian 
is composed of at least two terms: the vacuum contri- 
bution Hq and the MSW potential V. When solving 
for S one must work in a particular basis and the ba- 
sis determines the structure of the terms in the Hamil- 
tonian. In the 'mass' basis the vacuum Hamiltonian is 
diagonal and described by two mass squared differences 
Smfj = mj — m 2 and the neutrino energy E. Through 
this paper we shall use the values of Sm^i = 8 x 10 -5 eV 2 , 
|£m§ 2 | = 3 x 1(T 3 eV 2 , sin 2 28 12 = 0.83 and sin 2 26» 23 = 1 
which are consistent with present experimental values. 
In the flavor basis the off-diagonal elements are non- 
zero leading to the phenomenon of flavor oscillations. 
The two bases are related by the Maki-Nakagawa-Sakata- 
Pontecorvo [3!j|4(| unitary matrix parametrized by three 
mixing angles, 6*12, #13 and 6*23, a CP phase and two Ma- 
jor anna phases. 

In contrast, the MSW potential is diagonal in the fla- 
vor basis because the matter picks out the neutrino fla- 
vors. The common neutral current contribution to the 
MSW potential may be dropped because it leads only 
to a global phase which is unobservable leaving just the 
charged current potential, \^2G FTi e (r) where Gf is the 
Fermi constant and n e (r) the electron density, which af- 
fects just the electron neutrino/ antineutrino i.e. the el- 
ement V ee . In addition to the MSW potential, it has 
been found that the neutrino density in supernovae is 
so high that an additional potential due to neutrino 
self-interactions must be included. This neutrino self- 
coupling has been shown to lead to very interesting be- 
havior but for our purposes the self-interaction is negli- 
gible when the turbulent region in the star has moved 
beyond 1000 km so we shall ignore this contribution. 

When the vacuum and matter terms are added to- 
gether the Hamiltonian is neither diagonal in the mass 
nor the flavor bases so one would expect oscillations of 
both the flavor and mass probabilities. These oscillations 
are a source of potential confusion for any analysis. A ba- 
sis can be found which diagonalizes H for a given value 
of the electron density in the sense that there is a matrix 
U such that U^HU = K where K is the diagonal matrix 
of eigenvalues. This basis is known as the matter basis 
which becomes the mass basis (up to arbitrary phases) 
when the MSW potential disappears. The matter mix- 
ing matrix U which achieves this diagonalization depends 
upon the position through the star therefore dU / dx 
in general. The non-zero derivative of the matter mixing 
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FIG. 1: The MSW potential as a function of distance through 
a supernova taken from a hydrodynamical simulation. The 
vertical lines indicate the positions of the reverse and froward 
shock in the profile. The horizontal dashed-dotted line is the 
two-flavor resonance density for a 25 MeV neutrino with mix- 
ing angle sin 2 26 — 0.1 and mass splitting Sm 2 = 3 x 10 -3 eV 2 



matrix re-introduces off-diagonal elements into the mat- 
ter basis Hamiltonian which will lead to mixing between 
the matter basis states if they become large. We refer the 
reader to Kneller & McLaughlin [55] and Calais, Kneller 
& Volpe [HI for a more detailed description of the mat- 
ter mixing matrix. We shall report our results using the 
matter basis states throughout this paper. 

Next we must introduce the turbulent density pro- 
file though which the neutrinos will propagate. Ideally 
one would like to use density profiles taken from multi- 
dimensional simulations but at the present time that is 
not possible. The current multi-dimensional simulations 
do not extend out to the region of r > 10 4 km where 
the turbulence would have its greatest effects because 
the matter there has little bearing upon the explosion, 
and even if they did, they do not run to sufficiently late 
post-bounce times to see the shock move out there. Fi- 
nally, the dynamic scale the simulations would need to 
cover would be of order forty to fifty decibels - four to 
five orders of magnitude - because the neutrino oscilla- 
tion wavelength is significantly smaller than the radius in 
the high-density resonance region and beyond. For these 
reasons the effect of the turbulence upon the neutrinos 
is most often modelled as a random field. We adopt a 
one-dimensional supernova profile from a hydrodynami- 
cal simulation and in order to facilitate comparison the 
profile we select is taken from Kneller, McLaughlin & 
Brockman and is the same profile used in Kneller & 
Volpe. This profile is shown in figure ([1}. In the fig- 
ure we find two shocks: the forward shock at r s formed 
from the core bounce, and the reverse shock at ry formed 
by the wind created above the proto-neutron star run- 
ning into the material ahead of it. In multi-dimensional 
simulations of supernova both these shock fronts are dis- 
torted leading to strong turbulence in the region between 
them. But for this paper we shall use neutrino energies 
and mixing parameters such that the H resonance density 
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FIG. 2: The transition probability Pi 3 in an inverted hierar- 
chy (top panel) and P23 in a normal hierarchy (lower panel) 
after passing through the density profile shown in figure JT]) 
for a 25 MeV neutrino as a function of the mixing angle 813 



A. Modelling the turbulence 

The turbulence is introduced by multiplying the pro- 
file in the region between the reverse and forward shocks 
by a factor 1 + F(r) where F(r) is a Gaussian random 
field with zero mean. Since the quality of our results in 
this entire paper rests firmly upon us doing this well, it is 
worth our effort to explain carefully how F(r) was con- 
structed. The random field is represented using a Fourier 
series i.e. . 

F(r) = C* tanh (^—^J tanh (^^~ 

N k ' (2) 

x \/Vn {A„ cos (k n r) + B n sin (k n r) } . 

n=l 

for radii between r r < r < r s and zero outside this range. 
The two radii r r and r a are the positions of the reverse 
and forward shock respectively found in the underlying 
profile. In this equation the parameter C* sets the am- 
plitude of the fluctuations. The two tanh terms are in- 
cluded to suppress fluctuations close to the shocks and 
prevent discontinuities at r s and r r , and the parameter 
A is a scale over which the fluctuations reach their extent 
size. We set A = 100 km. In the second half of equa- 
tion |(2} the members of the set of co-efficients {A} and 
{B} are independent standard Gaussian random variates 
with zero mean thus ensuring the vanishing expectation 
value of F. The Nk wavenumbers form a set k n with 
power spectrum E{k) and, finally, the parameters V n are 
k-space volume co-efficients to which we return shortly. 
The power spectrum of the random field was selected to 
be 



does not intersect the shocks. The reason we avoid the 
shocks is twofold. Hydrodynamical simulations typically 
yield 'soft' shocks that do not cause transitions between 
the neutrino states if the mixing angle is too big. This 
lack of a transition is unphysical. The second reason is 
that we wish to focus solely upon the turbulence effect 
and diabatic MSW transitions caused by the shocks com- 
plicates the interpretation. For these reasons we will use 
25 MeV for the neutrino energy and the two-flavor reso- 
nance density for a 25 MeV is shown in the figure. The 
reader can verify that it does not intersect the profile at 
either shock. For future reference, the transition proba- 
bility P23 for a normal hierarchy and Pi 3 for an inverted 
hierarchy are very close to zero for this neutrino energy, 
and mass splitting of Sm 2 = 3 x 10 -3 and a mixing angle 
given by sin 2 2#i3 =0.1. The reader may be surprised to 
see that the figure indicates the transition probabilities 
change from the diabatic limit, P13 = 1 or P23 = 1, to the 
adiabatic limit P13 = or P23 = is not monotonic in 
the H resonance channel because of the aforementioned 
presence of the multiple H resonances in the profile. The 
multiple resonances leads to an interference effect which 
is sensitive to #13 when 0.1° < $13 < 1° for this neutrino 
energy, profile and mass splitting [Hj, [HJ . 



E(k) = 



(<*-!) 
2fc* 



(3) 



Here fc* is the cutoff scale, a is the spectral index and is 
the Heaviside step function. Throughtout this paper we 
shall use a wavenumber cutoff fc* set to fc* = n/ (r s — ry) 
i.e. a wavelength twice the distance between the shocks 
and we shall adopt the Kolmogorov spectrum where 
a = 5/3. The method of fixing the Nk k's, Vs, A's and 
S's for a realization of F is 'variant C of the Random- 
ization Method described in Kramer, Kurbanmuradov, 
& Sabelfeld 58]. This Randomization Method partitions 
the k-space into Nk regions and from each we select a 
random wavevector using the power-spectrum, E(k), as 
a probability distribution. The volume parameters V n are 
the integrals of the power spectrum over each partition if 
the power spectrum is normalized to unity. Variant C of 
the Randomization Method divides the k-space so that 
the number of partitions per decade is uniform over Nd 
decades starting from a cutoff scale fc*. The logarithmic 
distribution of the modes is designed to ensure the quality 
of the agreement between the exact statistical behavior 
of the field and that of an ensemble of realizations is uni- 
form over a the range of lengthscales considered i.e. it is 
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FIG. 3: The reduced wavelength of the splitting between the eigenvalues. The density profile is that shown in figure |T]), the 
neutrino energy is 25 MeV and the neutrino mixing angle #13 is given by sin 2 2#i3 = 0.1. In the left panel we show the normal 
hierarchy case, in the right panel the inverted and the top row of each is for the neutrinos and the bottom row for antineutrinos. 



scale invariant. This feature is important for our study 
because the oscillation wavelength of the neutrinos is con- 
stantly changing as the density evolves. The evolution 
of the reduced oscillation wavelengths for the neutrinos - 



A, 



l/\5ki. 



and antineutrinos - A,- 



l/\Sh 



where 



Skij and Skij are the differences between the eigenvalues 
i and j of the neutrinos and antineutrinos respectively - 
as a function of distance through the profile are shown 
in figure ([3]) for both a normal and an inverse hierarchy 
when the mixing angle #13 is set to sin 2 26*13 = 0.1 and the 
energy is E = 25 MeV. Again the reverse and forward 
shocks are indicated by the two vertical dashed lines. 
This figure can be used to determine a suitable value for 
iVrf because we observe that in the region between the 
shocks the typical wavelengths are > 1 km which is the 
minimum lengthscale we need to cover [5ll [59[. This 
is approximately 4 orders of magnitude smaller than the 



turbulence cut-off scale 1/fc* thus we deduce that we need 
to pick Nd > 4 to cover the necessary decades in k-space. 

With Nd determined we now seek a suitable value of Nk 
by requiring that the statistical properties of an ensem- 
ble of random field realizations closely match the exact 
properties for the field. The statistical property we com- 
pute is the second order structure function G2(Sr) given 
by 



G 2 (Sr) = (F{r + Sr) - F(r)) 2 



(4) 



where Sr is the separation between two radial points. The 
function G2(Sr) is related to the two-point correlation 
function B(8r) via G2(Sr)/2 = 1 — B(5r) and for the 
power spectrum we have adopted we can compute the 
two-point correlation function analytically to be 



B(5r) 



(a -I) 



(2 7 rfc^<5r) a " 1 {exp(^) r(l - a, 2tn ^ Sr) + exp r(l - a, -2m fc* fr)} 



r 



(5) 



where T(n,x) is the incomplete Gamma function. In 
figure ((3]) we show the ratio R(5r) of the numerically 
calculated structure function to the exact solution as a 
function of the scale k+dr when we use either Nk — 50 
wavenumbers spread over Nd = 5 decades or Nk = 90 
wavenumbers over Nd = 9 decades. The numerical cal- 



culation is the average of 30,000 realizations of the tur- 
bulence and the error bar on each point is the standard 
deviation of the sample mean. The figure indicates that 
the method we use to generate random field realizations 
reproduces the analytic results for the structure function 
very well and with high efficiency because good agree- 
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FIG. 4: The ratio R(Sr) of the numerically calculated struc- 
ture function to the analytic result as a function of k*Sr. The 
two curves in the figure correspond to {Nt, Nd} = {50, 5} (red 
solid) and {N k ,N d } = {90,9} (blue dashed). At every k+Sr 
we generated 30, 000 realization of the field and the error bar 
on each point is the standard deviation of the mean. 

merit between the statistics of the ensemble and the ex- 
act result requires just Nk/Nd = 10. In fact, like Kramer, 
Kurbanmuradov, & Sabelfeld |58j ] before us, we find even 
Nk/Nd ratios of just Nk/Nd ~ 2 — 3 are sufficient to give 
acceptable agreement. We re-assure the reader we shall 
stick with N k /Nd = 10. 

III. RESULTS 

With the construction of the random fields in place we 
can proceed to generate a turbulent profile and propagate 
neutrinos and antineutrinos through it. This construc- 
tion and propagation recipe is then repeated a minimum 
of one thousand times - sometimes much larger - to con- 
struct an ensemble of transition probabilities of size N . 
Once we have our sample we can then go ahead and com- 
pute means (Pij), variances cr^, etc. The hierarchy will 
be set to normal and we shall comment on how our re- 
sults translate to the inverted hierarchy. The neutrino 
energy will be fixed at E = 25 MeV, typical of super- 
nova neutrino energies. The turbulence effects - or lack 
of them - when using a value of 6*13 close to the present 
measurements was not fully explored in previous stud- 
ies so to make a connection with previous works, and to 
explain why a large value of #13 gives the results that 
it does, we shall consider multiple values of #13 in order 
to show what other possibilities would have produced in 
contrast. 



A. Small amplitude turbulence 

For small amplitude turbulence only the H-resonant 
channel is affected: mixing between states v-i and 1/3 for 
a normal hierarchy and states v\ and D3 for an inverted. 



Previously it has been found that effects could appear 
in the neutrinos even for turbulence amplitudes in the 
range 10~ 5 < C* < 0.1 when the mixing angle was set at 
sin 2 26*13 = 4 x 10~ 4 . If we allow the value of #13 to float 
then we find the normal hierarchy H resonance channel 
transition probability P23 can become more or less dia- 
batic. This can be explained from the behavior of the 
diabaticity parameter T23 [551 ] . which characterizes the 
degree of mixing in the H resonance channel for a normal 
hierarchy. This quantity is inversely proportional to the 
difference <5/c23 between the eigenvalues and k$ and 
proportional to the derivative of the matter mixing angle 
6*13. Increasing 6*13 increases the eigenvalue splitting and 
also makes the resonance 'wider' in the sense that the 
change between the limiting values of the matter mixing 
angle 6*13 occurs over a greater extent reducing the mat- 
ter angle derivative. Both effects decrease the diabaticity 
and, for these reasons, reaching the depolarization limit 
for P23 becomes more difficult if the domain of turbulence 
is fixed as is the case here. If the profile were changed so 
as to allow a larger turbulence region then eventually one 
should expect to reach the depolarization limit no matter 
what the mixing angle. A similar argument applies when 
#13 becomes small: now the diabaticity increases as #13 
decreases because the splitting between the eigenvalues 
at the resonance decreases and the transition occurs more 
rapidly. Either way, as #13 varies the distributions for the 
transition probability P23 will differ from the uniform dis- 
tributions seen in Kneller & Volpe leading to subsequent 
evolution of the expectation values and distribution vari- 
ances. This evolution with 8 13 is seen in figure §5§ where 
we plot the mean value (P23) and the standard deviation 
of the samples from a single emission point as a function 
of #13. The reader should compare the evolution of (-P23) 
in this figure with that in figure The mean value of 
(-P23) has an inflection region between 0.1° < #13 < 1°: 




e 13 n 

FIG. 5: The mean of the transition probability, (P23) (solid 
line) and the standard deviation 023 (dashed line) as function 
of the mixing angle #13. The turbulence amplitude is set to 
C* = 1% and {Nk,Nd} = {50,5} and the neutrino energy is 
25 MeV. 
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for the smaller values of #13 we see (P23) > 1/2, for the 
larger (P23) < 1/2 and for 13 ~ 9° the mean value of 
P23 is almost zero. One also observes how the sample 
standard deviation changes as #13 varies and sees that it 
is maximal at 023 = 0.28 for the range 0.1° < 6*i 3 < 1° 
and almost zero when #13 ~ 9°. This figure shows how 
the measurement of #13 has brought clarity to the issue 
of turbulence and supernova neutrinos. For #13 outside 
the range 0.1° < #13 < 1° the distribution of P23 is es- 
sentially a delta function at either zero or unity; for 13 
inside the range 0.1° < #13 < 1° the distribution is uni- 
form. Thus when #13 was unknown it was impossible to 
determine whether the effect of small amplitude turbu- 
lence was negligible or overwhelming. The measurement 
of a large value of #13 indicates it is the former and the 
result has consequences for the observability of spectral 
features in the next Galactic supernova burst signal. 



B. Large Amplitudes 



The neutrino mixing channels 
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FIG. 6: The frequency distribution of the transition proba- 
bility P31 for a turbulence amplitude is set to C* = 0.3. The 
neutrino energy is 25 MeV and we used Nk = 50, Nd = 5 
for the turbulence generator. The top panel used a value 
of 813 given by sin 2 20i3 = 0.1 and the bottom panel is for 
sin 2 20i3 = 4 x 10~ 4 . 



For large amplitudes, C\ > 0.1, the effects of turbu- 
lence are no longer restricted to the H resonance channel 
but appear in numerous places. The first effect worth 
noting is that the distribution of the H-resonance chan- 
nel transition probability, P23 in the case of a normal 
hierarchy, becomes independent of #13. This can be seen 
in figure (JTJ) where the reader will observe the evolution 
of the mean value of this transition probability (P23) as 
a function of C*. For the two values of #13 considered, 
the spread in (P23) at small amplitudes has disappeared 
by C* ~ 0.3. We also notice that around this same tur- 
bulence amplitude there begins the shift to three-flavor 
depolarization where (P23) = 1/3. 

In addition to the changes in the H-resonance chan- 
nel we also begin to observe mixing in the L-resonance 
channel, between vi and i>2 as the amplitude grows. 
This simultaneous mixing between v\ and 1/2 and Vi 
and V3 breaks HL factorization and Kneller & Volpe 
gave two examples were given that explicitly showed 
broken HL factorization. For the neutrino mixing pa- 
rameters we are using, the ratio of H and L resonance 
densities (using the two- flavor formula) is ph/pl = 
(£m| 3 cos2#i3)/(5mf 2 cos2#i 2 ) 90. This large ratio 
would seem to imply that we need fluctuations of order 
F ru 1 because only if F = —0.99 would the density 
fluctuation give ph(^ + F) ~ pl- Three effects soften 
this requirement: the L resonance has a large width - 
Apl/pl = tan 2^12 ~ 1, the density in the turbulence 
region can be much lower than the resonance density pu 
for the given neutrino energy - see figure (JTJ), and, finally, 
our choice of a Gaussian random field for the turbulence 
will ensure that large fluctuations will occur occasionally 
no matter what we use for the amplitude, larger ampli- 
tudes just make the extremal fluctuations more proba- 
ble. The clearest signature of broken HL factorization 



is a non-zero transition probability P31 because only if 
HL factorization is broken can we generate an effective 
mixing between v\ and v%. To see this we consider the 
S'-matrices for the case of factored HL resonances and 
broken factorization. The S-matrix for passing through 
one or several H resonances, Sh , has the genreal form 



1 

S H = I a H Ph 

-/?* a H 



(6) 



where an and fin are Cayley-Klein parameters. Simi- 
larly the S-matrix for L resonances, <Sl, is 



St 



a L p L 

-Pi < 
1 



(7) 



where an and fin are Cayley-Klein parameters for the L 
resonance. If HL factorization holds then the S-matrix 
which describes the evolution for the neutrino through 
the entire profile is S = SlSh- If all the L resonances 
occur after the H resonances then we find the transition 
probability P31 — \Szi\ 2 is identically zero. But if addi- 
tional H and L resonances occur, denoted by S' H and S' L , 
then the S-matrix describing the neutrino evolution is of 
the form S = S' L S' H S L S H and P31 = W h \ 2 \Pl\ 2 will be 
non-zero. In figure ^ we show frequency distributions 
of the transition probability P31 for neutrinos at two val- 
ues of #13 when C* = 0.3. The distributions are clearly 
non-zero for non-zero P 3 i as expected if HL factorization 
were broken. The distributions fall rapidly as something 
like inverse-power laws, /(P31) oc I/P31 with n ~ 2 or 
exponentials for this particular calculation. 

Other mixing channels which were previously delta- 
distributed for small turbulence amplitudes - such as P12 
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FIG. 7: The mean of the transition probabilities Pi 2 - top 
panel - Pi 3 - center panel - and P23 - bottom panel - through 
anisotropic turbulence as a function of C\ for antineutrinos 
emitted from a single point. Each curve corresponds to a 
different value of #13: sin 2 2(9i3 — 4 x 10~ 4 are squares joined 
by solid lines, sin 2 26>i3 = 0.1 are circles joined by dashed 
lines. 



and Pi 3 - also begin to possess similar inverse power- 
law/exponential distributions when C* > 0.1 and their 
means increase quadratically with C*. The evolution of 
these two transition probabilities is also shown in figure 
(J7J). The figure shows that (P12) at some fixed C* in- 
creases as #13 increases but (P13) decreases as #13 in- 
creases though, in both cases, the change is not very 
large. The anticorrelation between (P12) and (P13) is 
a reflection of the unitarity requirement that S^P^ = 1 
for a given i. 



2. The antineutrino mixing channels 

In addition to breaking HL factorization, large am- 
plitude turbulence induces effects in the non-resonance 
channels particularly v\ f-> vi regardless of the hierarchy, 
V\ <-> ^3 for a normal hierarchy and P23 for an inverted hi- 
erarchy. If we stick with considering the normal hierarchy 
case then we can compute the mean of the non-resonant 
transition probabilities P12 , P13 and P23 as a function of 
the mixing angle #13 and turbulence amplitude C* . These 
are shown in figure ©. Like (P12) and (P13) shown in 




0.1 0.2 0.3 0.4 0.5 



FIG. 8: The mean of the transition probability P12 - top 
panel - P13 - center panel - and P23 - bottom panel - as a 
function of C* for antineutrinos. Each curve corresponds to 
a different value of #13: sin 2 2#i3 = 4 x 10 -4 are squares and 
sin 2 2613 = 0.1 are circles. 



figure (O, the reader will observe that the three tran- 
sition probabilities grow rapidly with C+ reaching the 
levels of (P12) ~ (10%), (P13) ~ (1%) and (P 23 ) ~ (1%) 
at C+ ~ 0.5. Further comparison with figure (J7J) reveals 
(P12) ~ (P12) and the expectation values for the tran- 
sition probabilities P13 and P23 are both smaller than 
(P13) by roughly an order of magnitude and much more 
sensitive to $13. The expectation value for P13 varied 
by a factor of ~ 2 when #13 allowed to float, here P13 
and P23 change by an ^ 1 — 2 orders of magnitude when 
increasing from sin 2 26*13 = 4 x 10~ 4 to sin 2 26*13 = 0.1. 
This same sensitivity to 6*13 was explained in Kneller & 
Volpe as due to the proportionality of the antineutrino 
diabaticity parameter Ti 3 to the vacuum mixing angle. 
The current preference for #13 close to #13 ~ 9° indicates 
(P13) and (P23) can be of order a few percent if C* ~ 0.5. 



IV. SUMMARY AND CONCLUSIONS 

The effects of supernova turbulence upon the flavor 
composition of neutrinos that pass through it depend 
upon the numerous parameters that one needs to in- 
troduce. For a neutrino energy of 25 MeV and using 
a supernova density profile taken from a simulation 4.5 s 
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post-bounce, turbulence of amplitude C* = 1% only af- 
fects the H resonance mixing channel to any apprecia- 
ble degreee and then only for mixing angles in the range 
0.1° < 6*i3 < 1°. For presently favored value of $13 closer 
to 9° there is little effect of small amplitude turbulence. 
This result will remain valid for other similar neutrino 
energies which have their MSW resonances in the region 
where the turbulence is located. At later post-bounce 
times the turbulence region will move outwards and to 
lower densities affecting the neutrino L resonance. The 
phenomenology will be similar to the affect of turbulence 
upon the H resonance. The removal of turbulence effects 
upon the neutrinos at small amplitudes has important 
consequences for the prospect of observing signatures of 
collective and shock wave effects in supernova neutrino 
burst signals, which is explored in Lund & Kneller [60(. 

For the same neutrino energy and post-bounce epoch, 
the turbulence effects metastasize as the amplitude in- 
creases. In addition, the sensitivity to #13 in the H reso- 



nant mixing channel is lost and HL factorization becomes 
increasingly broken. For amplitudes of > 0.3, and a 
normal hierarchy the expectation values of the transition 
probabilities P12, P13, P23 and Pyi are of order 10% or 
greter; in an inverted hierarchy it is the transition prob- 
abilities P12, P12, P13 and P23 whose expectation values 
are of equivalent magnitudes. These channels are the 
most promising for observing the signatures of turbu- 
lence. 
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